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On the scent of the knee — air shower measurements with KASCADE 
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Detailed investigations of extensive air showers have been performed with the data measured by the KASCADE 
experiment. The results allow to evaluate hadronic interaction models, used in simulations to interpret air shower 
data. The all-particle spectrum of cosmic rays and their mass composition, as well as individual spectra for groups 
of elements have been reconstructed. The results suggest, the inee in the all-particle cosmic-ray energy spectrum 
is caused by a rigidity-dependent cut-off of individual element groups. 



1. INTRODUCTION 

The earth's atmosphere is continuously bom- 
barded by highly relativistic ionized particles, 
first discovered and named "cosmic rays" by 
V. Hess in 1912. Present-day experiments show 
the cosmic-ray energy spectrum extending up to 
more than 10^° eV. The flux spectrum follows a 
power law dN/dE oc E~'^ over many decades in 
energy. The most prominent feature is the knee in 
the spectrum around 3 PeV where the spectrum 
steepens from 7 « 2.7 to 7 ~ 3.1. The origin 
of cosmic rays is still under debate. Strong, rel- 
ativistic shock fronts expanding from supernova 
explosions are favoured by popular models for the 
acceleration of ionized particles. Such models ex- 
plain the particle acceleration up to energies of 
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about Z ■ 10^^ eV, with the nuclear charge Z of 
the particle. This upper limit coincides for pri- 
mary protons with the mentioned steepening of 
the spectrum, and the origin of the knee is re- 
lated to the upper limit of acceleration in several 
models. 

Since the charged particles are deflected in the 
interstellar magnetic fields, the only hint of their 
sources are their energy spectrum and the mass 
composition, or more preferable, the energy spec- 
tra of individual elements. Cosmic rays at en- 
ergies below 1 PeV have been directly observed 
by balloon-borne instruments at the top of the 
atmosphere or in outer space. At higher en- 
ergies, the steep falling flux spectrum requires 
large detection areas or long observation periods, 
presently only possible in ground-based installa- 
tions. These detector systems measure the sec- 
ondary particles produced by cosmic rays in the 
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atmosphere, the extensive air showers (EAS). 

To investigate the cosmic rays from several 
10"'^^ eV up to 10^^ eV the air shower experi- 
ment KASCADE ("Karlsruhe Shower Core and 
Array DEtector" ) ||l| has been built on-site at the 
Forschungszentrum Karlsruhe in Germany. The 
experiment detects the three main components of 
EAS simultaneously. A 200 x 200 m^ scintillator 
array measures the electromagnetic and muonic 
components. The 320 m^ central detector system 
combines a large hadron calorimeter Q with sev- 
eral muon detection systems . In addition, high 
energetic muons are measured by an underground 
muon tracking detector Q. 

2. INVESTIGATION OF HADRONIC 
INTERACTIONS 

In order to obtain the mass and energy of the 
shower-inducing particle, measured EAS observ- 
ables are compared with predictions of simula- 
tions. These simulations describe the develop- 
ment of an EAS in the atmosphere and, there- 
after, the signal response for all individual parti- 
cles hitting the detectors. A critical part in the 
simulation chain is the model used to describe 
the high-energy hadronic interactions, since the 
model has to extrapolate into kinematical and en- 
ergy regions not covered by present-day collider 
experiments. The Karlsruhe EAS simulation pro- 
gram CORSIKA 1^ provides several high-energy 
hadronic interaction models — HDPM, DPM- 
JET, NEXUS, QGSJET, SIBYLL, VENUS — 
based on different phenomenological descriptions. 
Below 80 GeV the codes GHEISHA and URQMD 
are available. One objective of the KASCADE ex- 
periment is to evaluate these models and to pro- 
vide criteria for their improvement. 

The hadron calorimeter is a valuable detector 
for testing the interaction models. The structure 
of the hadronic component is examined in en- 
ergy and spatial coordinates. Observables used 
include the number of hadrons as well as their 
energy sum, their lateral distribution and their 
energy spectrum, the energy of each individual 
hadron relative to the most energetic hadron in 
an EAS, the maximum hadron energy, and the 
spatial distribution of the hadrons. All observ- 
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Figure 1. Frequency distribution for a) the num- 
ber of reconstructed hadrons and b) their energy 
sum in EAS for two shower size intervals |q|. 



ables are investigated as functions of the number 
of electrons and muons as well as of the hadronic 
energy sum. 

An example of such investigations is given 
in figure |l|. Shown are frequency distribu- 
tions for the number of reconstructed hadrons 
above 50 GeV and the hadronic energy sum for 
two shower size intervals for EAS with shower 
core inside the calorimeter. KASCADE find- 
ings are compared with predictions of the mod- 
els QGSJET, DPMJET, and NEXUS for two ex- 
treme mass scenarios, the primaries being only 
protons or iron nuclei Q . Combining all observ- 
ables, it turned out that the model QGSJET de- 
livers the most reliable description This 
model is used for the analyses described in the 
following sections. 
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Figure 2. Mean logarithmic mass (In A) vs. muonic shower size Nf^ for different sets of observables 



3. ENERGY SPECTRA AND MASS 
COMPOSITION 

Detailed investigations have been performed 
with the KASCADE data to study the influence 
of different observables and hadronic interaction 
models on the inferred cosmic-ray mass compo- 
sition. 

An example is shown in figure ^. In a Bayesian 
analysis the mean logarithmic mass (In^) = 
^jTi^i, where is the relative abundance of 
elements with mass Ai, is calculated for differ- 
ent combinations of observables. Several electro- 
magnetic, muonic, and hadronic observables have 
been used, including number of electrons Ne and 
muons N^^, number of hadrons above 100 GeV 
N^^'^, their energy sum '^Eh, the maximum en- 
ergy of a hadron per shower the number 
of high energy muons (i?^ > 2 GeV) N*, and 
their spatial distribution D_q. The left graph 
summarizes results obtained including the elec- 
tromagnetic shower size iVe, whereas the results 
in the right panel are obtained without N^. The 
combinations including lead to smaller (In A) 
values as compared to those without. The to- 
tal spread between the results is A(lnA) « 0.8. 
Similar systematic effects have been observed ear- 
lier for different hadronic observables These 



discrepancies are most probably caused by incon- 
sistencies in the hadronic interaction models used 
to interpret the EAS data. 

Several analysis methods have been carried out 
to obtain the cosmic-ray energy spectrum. The 
KASCADE data show a knee in the electromag- 
netic, muonic, and hadronic shower size spectra 
||ic| ]. From the size spectra, the primary energy 
spectrum has been derived (see e.g. includ- 
ing the absolute flux spectrum obtained for the 
first time using hadronic observables fl^ . 

A recent analysis uses electromagnetic and 
muonic shower size spectra in three different 
zenith angle bins. With a four component as- 
sumption for the mass composition of primary 
cosmic rays (protons, helium, CNO-group, and 
iron group), an unfolding algorithm is applied, 
taking into account shower fluctuations and ex- 
perimental effects. The hadronic interactions 
have been simulated using the model QGSJET 
above 80 GeV and the GHEISHA code below. In- 
dividual energy spectra for the four mass groups 
are obtained as shown in figure I Each spec- 
trum shows a Jcnee-like structure. The energy 
dependence of these cut-offs suggests a rigidity- 
dependent behaviour. 

The mean logarithmic mass calculated from the 
individual energy spectra is shown in figure ^ in- 
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Figure 3. Cosmic ray energy spectrum for four 
groups of elements and the resulting all-particle 
spectrum [O]. 
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Figure 4. Mean logarithmic mass vs. energy cal- 
culated from individual spectra for four groups of 
elements as shown in figure || [|l3| . 



dicating an increase of the average mass above 
the knee. 

4. CONCLUSION 

The EAS observables obtained by KASCADE 
are sensitive to hadronic interaction models used 
in simulations to interpret EAS data. At present, 
the combination CORSIKA/QGSJET best de- 
scribes the measurements. The systematic de- 
pendence of the mean logarithmic mass on differ- 
ent observables and models has been investigated. 
The systematic error for the model QGSJET is in 
the order of A(ln A) « 0.8. The reconstruction of 
individual energy spectra for groups of elements 
indicate a rigidity dependent cut-off, which ex- 
plains the knee in the all particle cosmic-ray flux 
spectrum. 
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